Planar pyrochlore, quantum ice and sliding ice 
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We study quantum antiferromagnetism on the highly frustrated checkerboard lattice, also known 
as the square lattice with crossings. The quantum Heisenberg antiferromagnet on this lattice is 
of interest as a two-dimensional analog of the pyrochlore lattice magnet. By combining several 
approaches we conclude that this system is most likely ordered for all values of spin, S, with a Neel 
state for large S giving way to a two-fold degenerate valence-bond solid for smaller S. We show next 
that the Ising antiferromagnet with a weak four-spin exchange, equivalent to square ice with the 
leading quantum dynamics, exhibits long range "anti-ferroelectric" order. As a byproduct of this 
analysis we obtain, in the system of weakly coupled ice planes, a sliding phase with XY symmetry. 



PACS numbers: 75.10.Jm, 75.10.Hk 

Introduction: On the heels of recent progress in un- 
derstanding highly frustrated classical magnets, coupled 
with a substantial experimental effort, El renewed atten- 
tion is now focused on the behaviour of their quantum 
counterparts. In particular, the Heisenberg pyrochlore 
antiferromagnet is being studied with view to the ques- 
tion of whether frustration-enhanced quantum fluctua- 
tions might lead to unconventional ordering - or com- 
plete absence thereof - especially for small values of the 
quantum spin, £.l3i3 This model classically is special in 
that frustration prevents any sort of ordering or dynamk 
cal phase transition down to the lowest temperatures ,ETtil 
for which reason it is termed a cooperative paramagnet 
or classical spin liquid. 

The challenge of this problem arises from the small en- 
ergy scale generated by the frustration: in a semiclassical 
picture, any linear combination of the classically degen- 
erate ground states-^, the collection of which is of exten- 
sive dimensionality!!^ - may be selected as the quantum 
ground state. For the highly frustrated two-dimensional 
magnet on the related kagome lattice, exact diagonali- 
sations of small clustersE 2 ] have provided crucial bench- 
marks. This system has turned out to be particularly 
well suited to this approach as it appears to have a very 
short correlation length, although one does find a large 
number of low- lying singlet excitations. 

For the pyrochlore magnet, it looks as if such results 
will elude us for some time to come. The pyrochlore lat- 
tice, being three-dimensional, displays a more inclement 
scaling of the Hilbert space dimension with linear system 
size. Moreover, its unit cell contains four spins and its 
structure implies that the smallest system without spuri- 
ous boundary condition effects contains at least 16 sites. 

To evade this, attention has shifted to a system which 
avoids some of these complications, namely the checker- 
board lattice (Fig. [j]). It is expected to have similar prop- 
erties to the pyrochlore as it has the same local structure 
- both can be thought of as networks of corner-sharing 
tetrahedra. Further, the size and topology of its ground 
state manifold for Heisenberg magnets are identical to 



the pyrochlore case However, it has a unit cell of only 
two spins, is two dimensional, and has allowed exact di- 
agonalisations of a good number of finite size systems of 
up to 36 spins. Such diagonalisations hajze very recently 
been juarried out by Palmer and Chalkeit3 and by Fouet 
et al&3 and other workers have-alaa recently studied this 
system by several techniques. liEnlZru 




FIG. 1. Phase diagram of the checkerboard lattice Heisen- 
berg magnet as a function of spin, 5*. The plaquette (left - fat 
bonds are those with enhanced probability of singlet forma- 
tion) and Neel (antiferroelectric) phases. A projection of the 
3-d pyrchlore lattice gives the 2-d checkerboard lattice (mid- 
dle). A tetrahedron becomes a square with crossings upon 
projection, all bonds of which have equal strength. 

In this paper, we take a somewhat broader view of 
quantum antiferromagnetism on the checkerboard lat- 
tice. For the Heisenberg problem, apart from the small 
S cases, which we study by a dimer model analysis and 
an Sp(A^) mean field theory, we analyze the large S re- 
gion, both within Sp(A r ) - which is also able to treat 
the intermediate region - and through the semiclassi- 
cal 1/S expansion. Here we find a strong ordering ten- 
dency - into a two-fold degenerate valence bond crystal 
for small S and a Neel state at large S (Fig. Q). Our 
predictions find support in the numerical work of Fouet 
et al. as we discuss further below. The nature of the va- 
lence bond order is at odds with other recent work.cl The 
5 = 1/2 kagome magnetli 2 ! behaves very differently from 
the checkerboard - and therefore probably the pyrochlore 
- ones. Details of the ordering, however, depend on prop- 
erties of the checkerboard lattice (most prominently, the 
explicit breaking of the symmetry between bonds in a 
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tetrahedron and the existence of nontrivial closed loops 
of length four residing on a bipartite lattice) which it does 
not share with the pyrochlore lattice; there we expect the 
order to be much more delicate. 

In addition we consider the Ising antiferromagnet with 
a weak four spin exchange dynamics - the Ising ground 
states are isomorphic to those of square ice and the dy- 
namics represent the shortest ring exchanges in this man- 
ifold arising from quantum tunneling - hence "quantum 
ice". The resulting state displays long range "antiferro- 
electric" order in ice terminology, or Neel order in Ising 
terminology (Fig. |l|). A byproduct of this analysis is the 
state of a stack of weakly coupled ice layers, which is seen 
to beia a "sliding" phase of the kind discussed in recent 
workjia albeit one that is protected by much simpler ar- 
rangements. 

We turn now to the details of these assertions and 
sketches of the underlying analyses. We treat the Heisen- 
berg problem first. 

Sp(iV) generalization: The Sp(iV) technology for frus- 
trated magnets, .-a reformulation of Schwinger boson 
mean field theorytj controlled by the introduction oil/N 
as a small parameter, was introduced and described in 
detail by Read and Sachdev.E3 One starts by rewriting 
the SU(2)~Sp(l) spin operators in terms of bosonic op- 
erators, with the constraint of b].b^ + &J&j = 
2S on each site, and S z = (&}b T - &|&|)/2, S+ = 

The antiferromagnetic nearest-neighbour Heisen- 
berg Hamiltonian is rewritten in terms of the bosonic 
operators (up to a constant): H = JY^Uj) Si • Sj = 

(-J/2) E {l3) {^ T bl b }r)(^ u K b ^)- One generalises this 
expression to Sp(iV) by formally introducing N flavours 
of bosons, labelled by capital letters, and replacing the 
antisymmetric tensor in H by its Sp(iV) generalisation 
J~ab = e>iu ^AB- In the limit N — ► oo at a fixed boson 
number per flavour per site, re, one obtains to leading 
order in 1/N a mean-field theory for spin S — re/2; fluc- 
tuations about this give rise to a gauge theory. 

Self-consistent solutions to the mean field theory are 
obtained by minimising with respect to the Hubbard- 
Stratonovich link fields Qy = {JAgbf bfJ), subject to 
a constraint on the boson number re. In the process, 
one obtains a dispersion relation for the energies to of 
bosonic spin 1/2 quasiparticles. Within the mean- field 
theory, there is generically a disordered phase at small re; 
when one of the bosonic modes goes soft as re is increased, 
condensation of these 'spinons' occurs, and long-range 
spin order ensues. 

Our results at mean-field level are readily summarised. 
We find zero expectation value of the diagonal bond vari- 
ables at all re. On the remaining bonds, long-range Neel 
order (Fig. |l|) develops above re c = 0.393, exactly as for 
the simple square lattice. For re < re c , the Neel correla- 
tions are only short-ranged. 

Whether or not the spinon excitations in the disor- 
dered, small re phase remain deconfined can only be 



settled by going beyond mean-field theory. As on the 
square lattice, one obtains a compact U(l) gauge theory 
at 0(1 /N), in which instanton tunneling effects lead to 
the formation of bond (Peierls) order .Ell An analogous 
calculation is presented by Chung et al. for the Shastry- 
Sutherland model.EJ Details of the ordering pattern, in 
particular those due to the inequivalence of the plaque- 
ttes with and without crossings, are perhaps most easily 
studied through a quantum dimer model (QDM). 
Quantum dimer model: This approach starts from 
the assumption that the magnet is in a regime where the 
Neel state is destabilised and the effective degrees of free- 
dom are singlet bonds between neighbouring spins, also 
called valence bonds, which are represented by dimers. 
As each spin participates in exactly one singlet bond, 
the Hilbert space consists of all hardcore dimer coverings. 
An effective dimer Hamiltonian is obtained by means of 
an overlap expansion, described in Ref. It is for- 
mally perturbative in a small parameter arising from the 
non-orthogonality of the spin wavefunctions describing 
the different dimer coverings. To zeroth order, all of the 
exponentially numerous dimer coverings are degenerate, 
but at the next order, a resonance term is generated. 
One obtains a resonance loop of length 4 corresponding 
to flipping a pair of dimers by 90° (Fig. |): I I ^ ~. The 
question is what kind of quantum dimer state is selected 
by this resonance move. 

The most important difference between the well- 
studied square lattice QTJM and that on the checkerboard 
lattice is the followingo As the overlap expansion is es- 
sentially organised by length, L, of the resonance loops, 
the leading order terms arise from the shortest possible 
resonance loops, namely those of length four. Carrying 
out the leading order overlap expansion, we find that 
all the loops of length four on plaquettes with crossings 
(the projected tetrahedra) have zero kinetic energy, t, 
in the dimer model and hence play no role. By contrast, 
the resonance loops on plaquettes without crossings, have 
nonzero kinetic energy (Fig. ||). Note that such loops of 
length four are absent from the pyrochlore lattice, where 
the shortest loop contacting more than one tetrahedron 
has length six. j—. 

Extending the results from the square lattice QDM,E3 
one therefore expects the resulting state to be a valence- 
bond crystal; the staggering of t strongly favours the pla- 
quette crystal, with the dimers resonating on one of two 
sublattices of plaquettes without crossings (Fig. [l]) . Note 
that the degeneracy of this state is two rather than four, 
as would be the case for the plaquette state on the square 
lattice, because of the explicit symmetry breaking intro- 
duced by the presence of the plaquettes with crossing 
interactions. 

Semiclassics: We now consider the case of large spin 
S and large re for SU(2) and Sp(iV), respectively. In 
either case, one compares the zero-point energy due of 
the excitations (spin waves and spinons, respectively) of 
different, classically degenerate ground states. This, in 
principle, requires evaluating that energy for the entire 
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ground state manifold, which is of extensive dimension- 
ality. 

For Sp(iV), we have compared the zero-point energies 
of all four-sublattice states, as well as of eight-sublattice 
coplanar states, and found that the Neel state selected 
at k c is also favoured in this limit, which suggests its 
stability for all n > n c . 

For large S, we have computed the zero-point energy 
of all four-sublattice states. We find good qualitative 
agreement with Henley's suggestiorO of an effective en- 
ergy functional of the biquadratic type, — Yliij) (S< ' Sj) 2 , 
in that it correctly reproduces the location of the maxi- 
mum and the minima. There remains a clegeneracy be- 
tween some inequvivalent collinear states.a One of the 
remaining degenerate states is indeed the Neel state, but 
other states, disfavoured at 0(1/ k), still have exactly the 
same energy of zero-point fluctuations at 0(1/ S). 
Implications for the 3-d pyrochlore magnet: The 
actual ordering behaviour -of-the pyrochlore is still far 
from settled. At this stageoO'Bj it appears to be some- 
what closer to that of the checkerboard lattice than ei- 
ther of the two is to the kagome case, where there has 
so far been no strong indication of long-range order of 
any kind for S = 1/2, and where the large- S state is 
necessarily non-collinear. We emphasize, however, that 
details of the ordering we find, such as the pattern of 
the bond solid or the size of the unit cell, crucially de- 
pend on differences between the two lattices - such as 
spatial dimensionality, the presence of inequivalent links, 
and short closed loops linking different tetrahedra. More- 
over, we remark ^at-mpst approaches employed so far for 
d = 3 pyrochloreDD&Q explicitly remove the equivalence 
of all tetrahedra by weakening the bonds on half of them; 
from such a starting point, it would would seem rather 
difficult to restore this equivalence, which is necessary for 
obtaining the plaquette ordering we find. 
Finite-size diagonalisations: The agreement of the 
above calculations on a number of central points is reas- 
suring. Most importantly, they all predict ordered states 
with translational symmetry breaking and a two-fold de- 
generacy. For small S, they lead us to expect dimer order 
of the plaquette flavour, which gives way to a Neel state 
at intermediate S. How does this compare to exact diag- 
onalisations of S = 1/2 Heisenberg magnets? 

Palmer and Chalker, who have studied systems con- 
taining up to 24 spins, find no clearly identifiable degen- 
eracy. Rather, there appears to be a large number of low- 
lying singlet states with a small gap and a much larger 
gap to triplet excitations. They rule out Neel order but 
are inconclusive about translational symmetry breaking. 
Fouet et al. agree with their results, but have in addi- 
tion studied a system with 36 sites. There, they find 
a particularly low ground state energy, suggesting that 
the boundary conditions for this system size accomodate 
well the quantum ground state. For this system, there 
does appear a two-fold near degeneracy of the ground 
state, with the states being described by the wavevectors 



expected for our dimer crystal. 

Quantum ice and sliding ice: We turn next to the 
Ising problem. The ground states of the Ising antiferro- 
magnet on the checkerboard lattice require two up and 
two down spins on each tetrahedron (square with cross- 
ings) . The six such possible configurations on each tetra- 
hedron can be identified with-the allowed vertices of the 
six-vertex model as followsEj Divide the square lattice 
consisting of the crossings at the centres of the tetrahe- 
dra into sublattices A and B in the usual fashion. Orient 
the links coming out of sublattice A(B) inwards if the 
spin sitting on that link is up (down) and outwards if the 
spin in down(up) (Fig. |l|). As all vertices are weighted 
equally, the ground state manifold has the extensive en- 
tropy, (3/4)ln(4/3) per spin, of square ice.E3 We will 
shortly find it useful that the ice problem also has a 
height representation in which an integer valued height 
living on the dual lattice steps up (down) by one onxross- 
ing an in(out) arrow clockwise around any vertexHj 

As a matter of principle, the degeneracy of the ice man- 
ifold will be lifted by quantum effects which will involve 
(in six vertex language) closed loops of arrows that will 
reverse direction. Here we study the ordering produced 
by this dynamics. The simplest such process for square 
ice involves a loop around a single plaquette. Translated 
into the Ising spin representation on the checkerboard lat- 
tice it gives rise to a Hamiltonian, Hq, acting between 
ground states of the classical Ising model, 

H Q = -Tj2(a + a-(T+a + + h.c.) (0.1) 

p 

where p denotes a sum over non-crossed plaquettes of the 
lattice, and are the raising (lowering) operators of 

the Ising spins a z = ±1. To study this Hamiltonian it is 
convenient to use an imaginary time discrete representa- 
tion of the path integral via the standard Trotter-Suzuki 
procedure - in this case it yields a set of ferromagnet- 
ically stacked planes of the checkerboard lattice, with 
an implicit time continuum limit at large ferromagnetic 
coupling, K T J2n a n a n+ii between neighbouring sites in 
adjacent layers, labeled by n£3 

We digress briefly to study the imaginary time rep- 
resentation at weak coupling, K T <C 1, which could be 
interpreted as the classical statistical mechanics of a set 
of square ice planes with a potential interaction between 
the planes. In this limit the interaction competes with 
the entropy of the planes and we may ask if it prevails de- 
spite its weakness. To answer this we switch to the height 
representation description of the ice problem. Herein the 
coarse grained heights in a gi*en plane are weighted by 
the pseudo-Boltzmann factorJ23 

p[h(x)] ot e-^i f d2x{vh)2 (0.2) 

while the spins (arrows) are represented by V x h for 
their small momentum components and by e l7vh for mo- 
menta in the vicinity of (ir,ir). For uncoupled planes, 
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the distribution rinPl^nM] ^ s a nxe d point of a two- 
dimensional renormalization group (RG) transformation 
in the standard fashion. The coupling between the planes 
induces two perturbations in the height language: the 
first of these, <j>\ — J2 n I d 2 x(Vh n ) ■ (Vh n +i), is ex- 
actly marginal under the RG while the second <f>2 — 
J2 n I d 2 x cos[7r(/i„+i — h n )] is irrelevant as is readily ver- 
ified. It follows then that the system exhibits a sliding 
phase at weak coupling, whence sliding ice. 
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FIG. 2. Spin correlations on sublattice 1 in the x-direction 
from Monte Carlo on 48 x 48 x 16 plaquettes. The classical 
curve, K T — 0, and K T =0.1 coincide. Neel long-range 
order, albeit weak, obtains for K T = 0.2. Inset: Numbers 
label sublattices, arrow translation vectors. In the QDM, one 
type of resonance move has zero kinetic energy, t = 0. 

Returning to the strong coupling problem, it is clear 
that we need to look for a phase transition out of the slid- 
ing phase into a flat phase where the layers lock. To this 
end we have carried out a Monte Carlo simulation whose 
results confirm the existence of the sliding phase as well 
as of the flat phase (Fig. ||). From our RG analysis, this 
should happen via a second order transition driven by (f>2 
becoming relevant, although we have not explored this 
in detail. Translated back the the quantum ice problem 
our results indicate that the system exhibits long range 
order of the antiferroelectric (or F model) kind (Fig. [I]), 
in which it maximizes the density of flippable plaquettes. 

For the spins on the checkerboard lattice lattice this 
is Neel order. The quantum dynamics is the leading or- 
der effect of either a transverse field (in which-epntext the 
Neel phase was conjectured by us previously )E3 or an XY 
exchange of either sign. One may wonder whether this 
conflicts with our previous analysis in suggesting Neel or- 
der at the Heisenberg point. However it has been shown 
within spin wave theorytLZl that the Neel state on the 
checkerboard lattice is unstable for S < 1 so there is 
every reason to expect a phase transition en route. For 
such order by disorder phenomena this is by no means ex- 
ceptional; one generically finds that fluctuations generate 



an ordered state out of a disordered ensemble, which they 
then in turn destroy as their strength increases. 

In summary, we have explored quantum frustrated an- 
tiferromagnetism on the checkerboard lattice. Despite 
an enormous classical (S — oo) degeneracy, we find a 
robust ordering tendency for Heisenberg magnets of any 
spin; for stacked square ice (Ising spins), a sliding phase 
precedes an antiferroelectric ordering transition. 
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